Hydrodynamics of confined active fluids 
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We theoretically describe the dynamics of swimmer populations confined in thin liquid films. 
We first demonstrate that hydrodynamic interactions between confined swimmers only depend on 
their shape and are independent of their specific swimming mechanism. We also show that due to 
friction with the walls, confined swimmers do not reorient due to flow gradients but the flow field 
itself. We then quantify the consequences of these microscopic interaction rules on the large-scale 
hydrodynamics of isotropic populations. We investigate in details their stability and the resulting 
phase behavior, highlighting the differences with conventional active, three-dimensional suspensions. 
Two classes of polar swimmers are distinguished depending on their geometrical polarity. The first 
class gives rise to coherent directed motion at all scales whereas for the second class we predict the 
spontaneous formation of coherent clusters (swarms). 



Soft materials composed of motile particles have 
seen a surge of interest over the last couple of years. 
They encompass auto-phoretic colloids [1 , self-propelled 
droplets [2, and vibrated grains [3JH]- This interest was 
triggered by their fascinating structural and transport 
properties akin to the one found in biological systems 
such as bacterial suspensions, migrating cells, and cy- 
toskeletal extracts (see Ref. [5] and references therein). 
These so-called active fluids are ensembles of self-driven 
particles capable of propelling themselves in the absence 
of any external actuation [5H9]. From a theoretical per- 
spective, these systems are commonly separated into two 
classes depending on the way they exchange momentum 
with their surroundings [5H7]. "Dry" systems, typically 
walkers, or crawlers, achieve locomotion by transferring 
momentum to a rigid substrate, and interact via short 
range contact interactions. In contrast "wet" systems, 
typically suspensions of swimmers, conserve momentum, 
and the particles interact at finite distance via long-range 
hydrodynamic interactions. A number of experimentally 
relevant situations involve monolayers of active particles 
living in confined fluid films, and thus belong to both 
classes - e.g. bacteria swimming on the surface of a cell- 
culture gel, or active colloids and droplets moving in mi- 
crofluidic channels [SJ \TU\ [TT] . 

In this letter, we describe the phase behavior of ac- 
tive fluids confined in two-dimensional (2D) geometries. 
In order to do so, we first revisit the description of hy- 
drodynamic interactions under confinement. We demon- 
strate that the far-field flow induced by a swimmer does 
not depend on the specifics of its swimming mechanism. 
The notions of pushers and pullers for instance, preva- 
lent in three dimensions (3D), are not relevant in thin 
films |12[ 113] . In addition, on the basis of a prototy- 
pal microscopic model, we show that due to friction with 
the walls, confined polar swimmers are not only prone 



to align along the local elongation axis but with the flow 
field itself. We then exploit these new interactions rules 
in 2D to address the large-scale dynamics of confined pop- 
ulations of swimmers. We establish a novel set of hydro- 
dynamic equations for confined active films, which quali- 
tatively differ from the modified Leslie-Eriksen equations 
for active liquid crystals [5 . An investigation of the re- 
sulting phase behavior leads to the distinction between 
two classes of polar swimmers depending on their geo- 
metrical polarity. The first class (large-head), gives rise 
to the emergence of coherent particle motion along the 
same direction at all scales whereas for the second class 
(large-tail), we predict the spontaneous formation of co- 
herent clusters (swarms). 

Let us consider an ensemble of self-propelled particles 
confined in a thin film of a Newtonian liquid. We address 
strongly confined geometries where the particle height 
is comparable to the film thickness, h, see Fig. [I] (left). 
At scales larger than h, the fluid flow is characterized 
by the projection of the z-averaged velocity field in the 
[x, y) plane. Far from a swimmer, the projected flow field 
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FIG. 1: Left: Sketch of a confined suspension of active par- 
ticles swimming freely in the (x, y) plane. Right: Close-up 
on a single polar swimmer (see text for notation). The active 
particles are confined between two walls in the z-direction. 
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u(r, t) is potential 



u(r) = -GVn(r), 



(1) 



where II(r) is the pressure at r = (x,y). The Darcy 
factor G scales as G ~ /i 2 /?y |14j . 

How does confinement affect hydrodynamic interac- 
tions between swimmers? In unbounded fluids, the flow 
induced by a swimmer depends on the microscopic details 
of the propulsion mechanism |15tfI7] . In the far-field, this 
flow is often well approximated by a force-dipole singu- 
larity, with a ~ 1/r 2 spatial decay, and as such has been 
used in most theoretical models |13lH51IT5] . This descrip- 
tion results in the distinction between so-called push- 
ers (or extensile swimmers), and pullers (or contractile 
swimmers). They corresponds to /orce-dipoles having op- 
posite signs, and displaying different large scale dynam- 
ics |13 [ H5 1 H9]. When confined by solid walls, these flows 
are screened algebraically and decay as ~ 1/r 3 , while 
retaining their angular symmetry. This screening of hy- 
drodynamic interactions was shown to suppress generic 
instabilities which are the hallmark of isotropic pusher 
suspensions [5]. 

As it turns out, however, the two main consequences 
of confinement has actually been overlooked so far. Any 
multipolar stress distribution on the surface of the swim- 
mer actually yields only subdominant contributions to 
the flow in the far field. For any particle transport mech- 
anism (swimming, driving, advection) the far-field flow 
induced by a particle moving in a confined fluid has in- 
stead the symmetry of a potential scmrce-dipole and de- 
cays as ~ 1/r 2 |20H22| . The distinction between pushers 
and pullers is thus irrelevant under confinement. Irre- 
spective of the propulsion mechanism, the flow induced 
by a swimmer located at r = R(t) is defined by Eq. 
and by a modified incomprcssibility relation 



V-u(r) = -er-V<5(r-R(<)), 



(2) 



where the dipole strength is cr = a R(i) — u'°)(R(t)) 

where u^ ' is the velocity field in absence of the particle, 
and a scales as the square of the particle size (for a disk- 
shape particle, a is twice the particle area) [20 . The 
dipolar solution, u d (r|R(i), cr), of Eqs. (JTJ— ([2J is given, 
for a particle located at the origin, by 



u d (r|0,<r) 



1 



2tt r 



(2ff - I) • cr, 



(3) 



with f = r/|r| and I the identity tensor [SOJ 122 • This 
framework has proven to accurately describe the inter- 
actions between confined advected droplets even in con- 
centrated systems |2"T1 l2"3Tf2"5] . Importantly, the angular 
symmetry of u d is different from the one of a force dipole: 
it is a polar flow field displaying the same angular depen- 
dence as that of a force monopole under confinement [22 
despite the swimmers being self-driven. The reason for 



this apparent paradox lies in the continuous momentum 
exchange with the confining walls, via the shear flow in 
the thin films that lubricate the swimmer-wall contacts, 
see Fig. [T] 

The second important difference with 3D suspensions 
concerns hydrodynamic interactions between swimmers. 
In order to account for these interactions, we first estab- 
lish the equations of motion of an isolated swimmer in a 
arbitrary fluid flow. We focus on swimming bodies with 
polar shapes, as is the case for most motile cells. For 
a swimmer at position R(i) we denote p(t) its orienta- 
tion (|p| 2 = 1) and v s the magnitude of its swimming 
velocity along p. From symmetry considerations and at 
leading order in |Vu|, the equations of motion of a polar 
swimmer for {R(t),p(i)} take the generic form 

Ra = V s p a + H±(6 a p - p a P0)uf) + {p a pp)U[3, (4) 
Pa = ^ (<W ~ PaP/3)up + v' ' (6 a p - PaPfi) ( V-fUp )p 1 , (5) 

where fj,± (resp. /in) is a transverse (resp. longitudinal) 
mobility coefficient and v and v' are two rotational mo- 
bility coefficients. In unbounded fluids, we have v = 
and /i_L = = 1, and Eq. (|5| then corresponds to 
Jeffrey's equation commonly used to quantify the ori- 
entation of anisotropic particles with the flow-elongation 
axis |131 IT!)] . Conversely, confined suspensions offer the 
possibility of having a nonzero value for v. Instead of 
reorienting due to flow gradients, swimmers can reori- 
ent because of the flow itself, a new type of orientational 
dynamics which has not been considered so far. 

To provide insight into the conditions for nonzero val- 
ues of v, we derive the above equations of motion for 
a prototypal microscopic model (dumbbell). We show 
how the lubricated friction with the walls induce both 
anisotropic mobility (fi± fxii) and a direct coupling 
between the flow velocity and the particle orientation 
(y 7^ 0). Consider a rigid-dumbbell swimmer, com- 
posed of two disks of radius b\ (resp. 6 2 ) located at Ri 
(resp. R2), and connected by a frictionless rigid rod of 
length a {b\, 62} (see Fig. [T] right). The lubrication 
forces between a disk-shape particle and the solid walls 
hinder its advection by the fluid. Passive disks would 
be transported at a velocity Rj(t) = /Zjii(Rj) (i = 1, 2), 
where the mobility coefficient /ij is comprised between 
(fixed obstacle) and 1 (passive tracer). We also intro- 
duce the drag coefficients ac when a disk is pulled by an 
external force F in a quiescent fluid, it moves at a veloc- 
ity R;(£) = otiF. Let us now assume that the two disks 
would propel at a velocity vi^p when alone, and let us 
compute the swimming speed and mobility coefficients 
from Eqs. Q-(|5| for the dummbell. The displacement of 
each disk results from the competition between (i) self- 
propulsion, (ii) the advection by the external flow u(°\ 
(iii) the advection of the disk i by the dipolar perturba- 
tion induced by the motion of the disk j, u (Rt|Rj, <Xj), 
and (iv) the inextensibility constraint, R2 — Ri = ap. 
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At leading order in hi /a, these contributions yield the 
following equations of motion for the "head" (i = 2) and 
the "tail" (i = 1) of the swimmer: 

Rx = i4 >p + Mu (0 HRi)+u d (Ri|R2,^ 2 )]+aiT,(6) 
R 2 = ^°) p + M2 [ u (o) (R2)+u d (R2 | Rij(Ti)] _ a2Tj(7) 

where the tension T ensures the inextensibility condi- 
tion, p • (R.2 — Ri) = 0. Defining the center of drag of 
the swimmer as R = (aiR2 + 02Ri)/(oi + 02), Eqs. |6|- 
([7]) are readily recast into the form of Eqs. with a 

dumbbell velocity and mobility coefficients given at lead- 



ing order by v s 



,(0) 



+ 0((h/a) 2 ), Mx = Q 2 /ii(l-72) + 



ai/i 2 (l -71), M|| = o 2 /ii(l+7 2 ) +oiAt 2 (l + 71) and v = 
Ka*2 +M172) - (Mi +M27i)]/a, where y t = bf(^ - 1) /a 2 . 
We first see that the translational mobility coefficients, 
/f±,l| depend only on the anisotropy of the swimmer, and 
are independent of its geometrical polarity (they remain 
unchanged upon a 1 -H- 2 permutation). In addition, 
as /in < /ij_j a non-swimming dumbbell making a finite 
angle with a uniform flow field would drift at a finite an- 
gle from the flow direction. We also obtain that indeed 
v 7^ for polar swimmers. Since the /ij's are decreasing 
functions of the particle radius, v is negative for large- 
head swimmers (62 > bi). and positive otherwise. From 
Eq. (5) we thus get that in a uniform flow large-head 
swimmers would reorient against the flow and propel up- 
stream. In contrast, large-tail swimmers (b\ > 62) would 
swim downstream. For apolar swimmers, v vanishes and 
the orientation of a symmetric dumbbell evolves accord- 
ing to the Jeffrey's orbits, Eq. ([5]), where v = and 
v' = a[fi 2 (1 +7i) + Mi(l + 72)]/2. Note that since u is ir- 
rotational, the orientation of an isotropic swimmer made 
of a single disk is not coupled to the background flow. In 
the rest of the paper we discard the conventional v' con- 
tribution to the orientational dynamics. It only yields 
short-wavelength corrections to the large-scale descrip- 
tion of polar-swimmers suspensions described below. 

We now turn to the dynamics of a dilute population of 
interacting swimmers in a quiescent fluid. We introduce 
the one-point probability distribution function, *&(]*, p, t) 
for swimmers with orientation p at position r and time t. 
The dynamics of the active particles is defined by Eqs Q- 
([5]), with the fluid velocity field, u(r, t), resulting from 
the linear superposition of force dipoles induced by each 
swimmer, u(r, t) = Jdpdr' ^(r', p, i)u d (r|r', er'), where 
a' = av s p. Assuming, that swimmers are subject to 
translational and rotational diffusion, ^(r, p, t) obeys the 
continuity equation 

dt% = -V • (*R) - V p • (*p) + L>V 2 * + A-Vp*, (8) 

where R and p are defined by Eqs. Q-Q, D and D r 
are the translational and the rotational diffusion coeffi- 
cients respectively, and V p stands for the gradient on 
the unitary circle. For simplicity, we neglect the trans- 



lational diffusion. Specifically, anticipating on our re- 
sults, we assume D <C v%/Dr, which is true for most 
biological and artificial micro-size swimmers. Note that 
for homogeneous suspensions, and due to the symme- 
try of the dipolar coupling, the sum of all hydrodynamic 
interactions vanishes: when V^r, p, t) = 0, we have 
Jdr'u d (r|r',cr') = 0, and thus from Eqs. Q-© it fol- 
lows that p = 0, and V • R = 0. The dynamics of an 
homogeneous population, from Eq. ([8]), reduces thus to 
the orientational diffusion of an isolated swimmer, and 
homogeneous phases relax toward an isotropic state over 
a time ~ . 

We now investigate the dynamic response of the ho- 
mogeneous and isotropic phase to spatial fluctuations 
of the concentration and orientation of the active par- 
ticles. The phase behavior is described in term of (i) 
the concentration field, e(r, t) = /^(r, p, t)dp, (ii) the 
local polarization, P(r, t) = - fpty(r, p, i)dp, and (iii) 
the local nematic-orientation tensor, Q(r, t) = - J(pp — 
gl)^'(r, p, t)dp. To establish their equation of motion, 
we need to add a closure relation to Eq. (|8| . As we focus 
on deviations from isotropic and homogeneous states, we 
expand ^ linearly in its three first moments [181 [T5] 

*(x, p, t) = — c (1 + 2 Pa P a + Ap a ppQ a p) , (9) 

where the numerical coefficients are chosen so that c, P, 
and Q are defined in a self-consistent fashion. Defin- 
ing \x = 1(//|| + fi±), and ft, = (/in — /tj_), and after some 
elementary but tedious algebra, the three nonlinear equa- 
tions of motion are inferred from Eqs. ([8])-([9]) as 



d t C = -V a [v s cP a + lLCU a + ftcQ aj 3Up\ 



(10) 



d t {cP a ) = 2 U » C ~ vc u pQf3a - D R cP a - VpZp a , (11) 
dt{cQ a p) = ^cu 7 (2<5 7 ( Q F ) g) - 6 a pP 7 ) - 4:D R cQ a/3 - V 7l 7 7Q ^, 



where the (potential) fluid velocity satisfies 



d a u a 



-<rv s d a (cP a ) , 



(12) 



(13) 



and where the expressions for the fluxes X and are 
given in supplementary information. 

Equations ( 10 )-( 13 ) fully describe the dynamics of 
the isotropic phase. We investigate their linear sta- 
bility with respect to plane-wave excitations of the 
form (5c, 5P, 6Q) cxp(ik • r — iu>t), with k = fcx. 
At linear order, we can integrate Eq. (13 1 for the 
fluid velocity, and recast the equations of motion 
into a set of two uncoupled linear systems hav- 
ing the form d t (5Py,5Q xy ) = M bcnd (6P y ,6Q xy ) and 
d t (Sc, 6P X ,5Q XX ) = M sp i ay (Sc, 5P X , 5Q XX ). The first sys- 
tem couples the transverse-polarization and the bend 
modes only. These modes are stable for all k, they cor- 
respond to damped sound-waves. The associated disper- 
sion relation is deduced from the eigenvalues of Afb cn d 
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FIG. 2: Stability diagram of a nearly isotropic and homoge- 
neous population of polar swimmers; Pe < (resp. Pe > 0) 
refers to large-head swimmers (resp. large-tail swimmers). 



as iuj — h(5Dn ± i\J —QD 2 R + (kv s /2) 2 ). In contrast, 
long-range hydrodynamic interactions between swimmers 
can destabilize the concentration (c), the longitudinal 
polarization (P x ) and the splay modes (Q xx ). To con- 
vey an intuitive description of this instability we intro- 
duce the two governing dimensionless numbers. First, 
Pe = vco(tv s /(2Dh) is a Peclet number comparing the 
rotational-diffusion rate Dr to the rate of rotation of 
a polar swimmer induced by a source dipole of magni- 
tude aC()V s (co being to the average concentration); large- 
tail swimmers (resp. large-head swimmers) correspond to 
Pe > (resp. Pe < 0). The second dimensionless num- 
ber, H = (flacQV s )/v s , compares the swimming speed, v s , 
to the advection velocity induced by a source dipole of 
magnitude <jcqv s . In the long- wave-length limit (k — ► 0), 
the eigenfrequencies associated with the stability matrix 
A^spiay then take the form 



oj c = 

LUp = 
WO = 



k 2 



l-H 
2D R Vl + Pe, 
-iD R (l + Pe) + 0(k 2 ), 
-AiD R + (D{k 2 ). 



(14) 

(15) 
(16) 



At th order in k, the total number of swimmers be- 
ing a conserved quantity we have w c = 0, and M sp i ay 
has only two non-trivial eigenvalues. Whereas rota- 
tional diffusion always stabilizes the nematic orientation 
{—iuiQ < 0), hydrodynamic interactions can in fact desta- 
bilize the isotropic state. From Eq. (15), we see that 
large-head swimmers with Pe < — 1 experience a generic 
instability: fluctuations of the local polarization are am- 
plified when the rotation induced by the hydrodynamic 
couplings overcome the diffusional relaxation of P x (see 
Fig. [2). 

Several comments are in order. First, although the 
growth rate of the instability does not dependent on k, 
the total polarization (k = 0) is not unstable. As dis- 
cussed above, the sum of all the hydrodynamic interac- 
tions cancels in this limit and no global directed flow 
can emerge spontaneously from an isotropic suspension. 



The instability shows however that groups of particles 
swimming coherently along the same direction form at 
all scales. Second, the generic nature of the instability 
is specific to the dipolar symmetry of the hydrodynamic 
interactions, and the polar shape of the particles, and 
can be intuitively rationalized as follows. From Eq. (13 1 



we see that any finite wave-length perturbation of P x 
along x results in a fluid flow in the opposite direction, 
with amplitude ~ ucqv s 8P x . Polar swimmers align with, 
or against, the local flow direction depending on their 
polarity. Large-head swimmers align along — u, thereby 
increasing the initial perturbation of P and destabilizing 
the isotropic state. Conversely, large-tail swimmers align 
in the opposite direction and the local polarization re- 
laxes to zero. As the reorientation rate of the swimmers 
is set by the magnitude of the velocity only (and not by 
the local strain-rate tensor), the growth (or relaxation) 
rate of the polarization is independent of the wave vector. 

This novel generic instability is qualitatively different 
from the one observed in unbounded suspensions of push- 
ers which, in contrast, is suppressed by confinement [S]. 
They differ in both the physical mechanisms at work and 
the structure of the unstable modes (bend versus splay 
modes). The only similarity is that in both systems the 
generic instability is a genuine collective effect due to the 
long-range nature of hydrodynamic interactions. 

To investigate the stability of the active film when 
Pe > — 1, we need to consider the eigenfrequencies , and 
the eigenmodes of M sp i ay up to 0(k 2 ). From Eq. (14 1 



we see that the combination of self-propulsion and rota- 
tional diffusion yields an effective diffusive dynamics of 
the suspension scaling as uj c ~ (v 2 / D R )k 2 , as could have 
been anticipated from the single swimmer problem |26| . 
However, hydrodynamic interactions result in a renor- 
malization of this single-swimmer effect. These inter- 
actions control both the magnitude and the sign of the 
effective translational diffusion. In the regions (Pe > — 1, 
H > 1) and (Pe < —1, H < 1), the effective diffusivity is 
negative and thus slowly destabilizes the isotropic phase 
(Fig. [2j. The associated eigenmodes are now complex 
superpositions of c, P x . and Q xx , and thus clusters of 
aligned particles form and propel in a coherent fashion 
(swarms), from a homogeneous film. Notably, both large- 
head (—1 < Pe < 0) and large-tail (Pe > 0) swimmers 
are prone to this second splay-destabilization mechanism. 
In the other regions of Fig. [2j the effective diffusivity is 
positive and concentration fluctuations are stable. 

In summary we revisited the theoretical description 
of confined populations of micro-swimmers. We showed 
that active particles interact hydrodynamically in generic 
manner, which is independent of the microscopic details 
of their propulsion mechanism and that, depending on 
their polarity: they may reorient in flows instead of solely 
flow gradients. Focusing on polar swimmers, we then 
constructed a large scale hydrodynamic theory from a 
minimal microscopic model (dummbclls). Our analysis 
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showed that the macroscopic orientational dynamics is 
very different from the modified Leslie-Eriksen model of 
active liquid crystals due to a difference in the symme- 
try of the microscopic coupling between confined polar 
particles and the fluid flow. It results in a novel phase 
behavior for active films and, in particular, spontaneous 
large-scale directed motion and swarming can emerge out 
of isotropic populations of confined swimmers. 
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Supplementary informations: Using Ricci-calculus notations, the expression of the fluxes X and J are: 



a 



1 

4 



flcP 7 u 7 S al 3 + (Ap, - jl)cP a Ufj + jlcPfjUa + Avcv s {Q afj + -^-) 



-if-Pj] - ^[Guydafj - U( a 5 i8 ) 7 ] + 2fj,u 7 cQ afj 



(17) 



(18) 



